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Abstract 

We study the simple gauge invariant model f 2 FF as a way to gen¬ 
erate primordial magnetic fields (PMF) in Natural Inflation (NI). We 
compute both magnetic and electric spectra generated by the f 2 FF 
model in NI for different values of model parameters and find that 
both de Sitter and power law expansion lead to the same results at 
sufficiently large number of e-foldings. We also find that the necessary 
scale invariance property of the PMF cannot be obtained in NI in first 
order of slow roll limits under the constraint of inflationary potential, 

V (0) 0. Furthermore, if this constraint is relaxed to achieve scale in¬ 

variance, then the model suffers from the backreaction problem for the 
co-moving wave number, k < 8.0 x 10 _7 Mpc _1 and Hubble parameter, 

Hi > 1.25 x 10 -3 Mpi. The former can be considered as a lower bound 
of k and the later as an upper bound of Hi for a model which is free 
from the backreaction problem. Further, we show that there is a nar¬ 
row range of the height of the potential A around A m i n ps 0.00874Mpi 
and of k around k m [ n ~ 0.0173Mpc _1 , at which the energy of the elec¬ 
tric field can fall below the energy of the magnetic field. The range 
of k lies within some observable scales. However, the relatively short 
range of k presents a challenge to the viability of this model. 
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1 Introduction 

Inflationary cosmology has solved many of the fundamental problems of the 
Big Bang model [1], Inflationary theory provides an explanation for the large 
scale structure of the Universe, which is linked to the quantum fluctuation of 
the field of inflation <p. The nearly flat, homogeneous, and isotropic universe, 
and the scale invariant spectrum of the anisotropies of the CMB, as well as 
the generation of primordial gravitational waves (PGW), are all predictions 
of inflation. All of the above predictions (with the possible exception of 
PGW) have been confirmed by various cosmological observations. 

Similarly, magnetic fields are being observed in all kinds of galaxies and 
cluster of galaxies at wide range of redshifts. Moreover, a lower bound, 
B > 3 x 10 -16 G, has been reported for intergalactic magnetic fields [2], 
Planck, 2015 results constrains the upper limit to be of the order of B < 
10 _9 G [2]. The commonly accepted model for the generation of galactic 
magnetic fields is the galactic dynamo [3]. This model needs a seed field to 
start amplifying the magnetic field. Also, it cannot explain the existence of 
fields in the absence of a uniform rotating charged medium, like the magnetic 
fields indicated in voids. Hence, such a detection may indicate a cosmological 
source of primordial origin. The generation of such a large scale primordial 
magnetic fields (PMF) is still an open question, see the reviews of this subject 
[4], However, one of the most interesting cosmological models is the simple 
f 2 FF model [5]-[7]. It has gained more interest because it is a stable model 
under perturbation. (See [8] and references therein.) Also, it can lead to a 
scale invariant spectrum of PMF under the exponential potential model of 
inflation [7]. 

The main difficulties with this model are the backreaction problem, where 
the scale of the associated electric field with PMF spectrum can exceed the 
scale of inflation itself [7, 9], and the strong coupling between electromagnetic 
fields and charged matter at the beginning of inflation [10]. In this case, 
if the electromagnetic field couples to charged matter, the physical charge 
associated with it is huge at the onset of inflation. For example, for the 
number of e-folds of inflation, N = 60, the physical electric charges, q oc e 120 
[10]- Such an absurdly large charge makes this model not viable. 

The starting point is the Lagrangian of a scalar (inflaton) field (t> coupled 
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to the electromagnetic (vector) field A /t [7] which can be written as, 


£ = -V-9 


\ (drf) (d^) + V (0) + (0, t) F m F vfi 


( 1 ) 


where, F u p = d„Ap — dpA u is the electromagnetic field tensor and g is the 
determinant of the spacetime metric g^ u . The first term in the Lagrangian 
is the standard kinetic part of the scalar field, and the second term,]/ ( 0 ), 
is the potential. A Lagrangian of a pure electromagnetic field would be of 
the form — jF^F^, but here we couple it to the scalar field through the 
unspecified function /(0, t). The main reason behind this coupling is to 
break the conformal symmetry of electromagnetism and hence prevent the 
dilution of the seed of magnetic field as it is generated in the inflation era. 

Simple models of inflation potential are based on a single scalar field, 
such as quadratic, V ( 0 ) ~ 0 2 , quartic, V (0) 0 4 [11, 13], Higgs potential, 

v(<t>) M 4 ^1 — exp [— [ 12 ] and the exponential potential, 

V (0) ~ exp[—y / 2 ii (0 — 0o)] [13]. The last one is used in [ 6 , 7] to find the 
magnetic and electric spectrum in the f 2 FF model. These models became 
more interesting after WMAP9 [14] and Planck [15]. As a result, the preferred 
potential class is the so called "plateau inflation ", at which V ( 0 ) 7 ^ 0 . 

In March 2014, results from the Background Imaging of Cosmic Ex- 
tragalactic Polarization (BICEP2) experiment were released [16]. They re¬ 
ported the detection of the tensor mode (B-mode) polarization of tempera¬ 
ture anisotropies in the CMB. The tensor to scalar ratio reported by BICEP 2 
was r = O. 20 q;o 5 , with r = 0 disfavored at 7.0<r. Also, the scale of inflation 
energy is close to the Grand Unified Theory GUT scale, Pgut 1 ^ 4 10 16 GeV. 
The results of BICEP2 put many standard inflationary models in trouble 
[17]- As a result, non-standard models, such as Large Field Inflation (LFI) 
[11, 18], Natural Inflation (NI) [19] models are a better fit with the BICEP2 
results. 

On the other hand, Planck results released in Sep 2014 [20] cast serious 
doubts on the primordial origin of this polarization signal. These results 
of Planck indicate that there is a significant contamination by dust over 
most of the high Galactic latitude sky in the same region where BICEP 2 
detected B-mode polarization. Consequently, there is a good chance that the 
source of the observations reported by BICEP 2 is all galactic dust and not 
of primordial origin. 

Finally, the joint analysis of BICEP2/Keck Array and Planck (BK/P) 
data was released on Feb 2015 [21]. The joint data of three probes eliminate 
the effect of dust contamination and show that the upper limit of tensor 
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to scalar ratio, ro .05 < 0.12 at 95% CL, and the gravitational lensing 13- 
modes (not the primordial tensor) are detected in 7a[21]. Similarly, the 
inflationary models are constrained based on the new analysis of data. The 
scalar spectral index was constrained by Planck, 2015 to be n s = 0.9682 ± 
0.0062(22], Consequently, the more standard inflationary models, like R 2 - 
inflation, which result low value of r, are the most favored one by Planck, 
2015. However, the chaotic inflationary models like LFI and NI are disfavored 
[22]. These results ruled out the first results of BICEP2, 2014. 

In this paper, the simple inflation model, f 2 FF, of PMF will be investi¬ 
gated in detail under the natural inflation (NI), in the same way as done in [7] 
and as we did in the context of LFI [23] and 7? 2 -inflatioii [24], We adopt nat¬ 
ural units, [c = h = ks = 1], the signature (—1,1,1,1), and a spatially-flat 
universe, where the reduced Planck mass, M P \ = (8tt<S') — 1 /2 ,will be taken as 
unity in the computations. Hence, the potential of NI can be written [12, 19] 
as, 

r«,)=A->[l + cos(b], ( 2 ) 

(J 

where A is fixed by CMB normalization of the potential, and a is the mass 
scale of the model. If a is taken to be of the order of M P i, then A 10- 13 M P1 . 
Since Eq.(2) is a periodic, even function of <fi, one can study the potential 
in the interval 0 G [0, na] [12]. In order to reach the GUT scale of inflation 
(10 15 — 10 16 GeV), the mass scale has to be of the order of a ~ M P \ [19]. 

The order of the paper will be as follows. In section 2, the slow roll 
natural inflation formulation will be presented for both the simple de Sitter 
model of expansion and the more general, power law expansion. In section 
3, the PMF and associated electric fields are computed for NI at different 
values of the parameters. In section 4, we summarize and discuss the results. 

2 Slow roll analysis of Natural Inflation 

During inflation, we will assume the electromagnetic field to be negligible 
compared to the scalar field, [7]. Hence, the equation of motion derived 
from (1) for the scalar field can be written as, 

(p + 3 Hcj) + 1% = 0, (3) 

where, H (t) = a ( t ) /a(t), is the Hubble parameter as a function of cosmic 
time, t, and a (t) is the cosmological scale factor. The over dot indicates 
differentiation respect to cosmic time, and 1% = 0 ( ,y. The Friedman equation 
can be obtained from the Einstein field equations by assuming a spatially-flat 
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Friedmann-Robertson-Walker FRW universe, which yields, 


H 2 


1 

3Mpi 2 


\$ + V W 


( 4 ) 


During inflation, the slow roll approximation allows us to neglect the second 
derivative in (3), which leads to the attractor condition 


0 


Y*. 

3 H' 


(5) 


Defining the slow roll parameters of inflation in terms of the potential 
[12, 25], of a single inflation field for NI (2), 


eiv ((/>) 



1 ( sin (*) V 
2 C 2 y 1 + cos(^) J 


e 2v (0) — 2Mpi 2 




2 1 

C 2 cos(^) + l’ 


( 6 ) 

(7) 


where ( = a/ Mpi. These slow roll parameters can be written in terms of the 
Hubble parameter, 


ei h (0) — 2M] 


pi 


Ho 

H 


, ^2 H ((f>) — 2M] 


PI 


H, 


<t>4> 


H 


The relation between the two formalisms [25] can be written as 

2 


Du — D h 


3 — C2 H 
3 — 6ih 


( 8 ) 

(9) 


All of the above parameters are assumed to be very small during the slow 
roll inflation,(eiy, e 2 vTmT 2 ir) W 1. Further, inflation ends when the values 
of {eip,ei h} ~> 1- At first order one can neglect and 6 2 f/as compared 
with 3, obtaining eiy ~ e\n- Therefore, using ( 8 ) and the relation between 
the cosmic time, t, and the conformal time, 77 , dt = a (r/) dr/, one can write 
the relation between conformal time and slow roll parameter, e\ h, as [18], 


V 


1 

aH + 


tlH 

~a?H 


da. 


( 10 ) 


Assuming, f\n ~ const, and then integrating (10) yields the power law ex¬ 
pansion of the universe during inflation, 


a{v) = l o\v\ 1 eiH , 


( 11 ) 
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where, Iq is the integration constant. 

On the other hand, in the simplest form of inflationary expansion (de 
Sitter), the universe expands exponentially during inflation at a very high 
but constant rate, Hi, 


Hi = - ~ const, 
a 


( 12 ) 


a(t) = a (■ t \) exp [Hit ], (13) 

where t\ is the starting time of inflation. In conformal time, Eq.(13) can be 
written as, 

a ^-~Hhr (14) 

Plugging (14) into the relation between cosmic and conformal time and in¬ 
tegrating implies that r; —» (—oo, 0 _ ) as t —> (0, oo). 

As inflation ends when e\y ((f)) ~ 1, then from Eq.(6), 

(pend ~ uarccos ( l + 2^2 ) • ( 15 ) 

If (pend <C 1, then a end <C 1, or a en d ^ (l - (2w 2 +1) 7rj, where, n =integer. 
Also, the relations between the slow roll parameters and the scalar spectral 
index, n s , and tensor-to-scalar ratio, r, can be written as follows [25], 


Tl s — 1 — 6 eiy + 2e2 V, 

(16) 

r = 16ei V . 

(17) 

Substitution of (6)-(7) into (16)-(17) yields 


3 1 2 f <P\ 

n,= l + - 2 -- 2 s e c 

(18) 

8 ( sin(-) \ 

(19) 

r "C 2 Vl4*cos(f)J • 


One can find the relation between r and n s which depends on the number 
of e-folds of inflation, N. The first order of approximation for N can be 
written [12, 25] as, 


N « 




1 

2Mpi 2 



v^i 


C 2 In 


1 - cos(^) 
l-cos(l) 


( 20 ) 
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where N is the difference between the final e-fold and the e-fold at t. Solving 
for p in (20) and substituting p en d from (15), 


p = aaiccos 


[1 — cos( 


'Pend 

G 


)] 


exp(-lV/C 



Now, combine (18) and (19), 



n q 


1 + sec 2 

c 



( Mt) V 

\{l + cos(£)}J ' 


( 21 ) 


( 22 ) 


By substituting (21) into (22), one can plot r — n s for different values of £ 
and N. 

Shortly after the onset of inflation the value of H becomes very high and 
is approximately constant, but later on it decreases as the value of the field 
changes. Also, after few iV*, the spacetime (pivot scale, &*) exits from the 
Hubble horizon. We adopt Planck, 2015 pivot scale, k* = 0.05Mpc _1 . One 
can consider H as a constant, after the first few e-foldings. That is basically 
the de Sitter expansion, which is exactly exponential expansion as described 
by (13). So, it is worthwhile to investigate both cases; the de Sitter, and the 
more realistic power law model described by (11). 


2.1 NI under de Sitter expansion 


The de Sitter model is the zeroth order approximation, which does not have 
graceful exit from inflation [26]. But it can be used as an approximation at 
the early stages of inflation. In conformal time, ?/, Eq.(5) can be written as, 


1 

a(v) 


V*_ 

3 Hi’ 


(23) 


where, Substitution of (2) and (14) into (26) and then integrating 

both sides yields, 

ac8c(-)d<f>= —. (24) 

V V 3 v 

Solving for p (rj), 

A 4 

p (rj) = 2cr arctan[c 2 ?y 3H i 2<r2 ], (25) 


where C 2 is the integration constant. By using the model of inflation at which 
the potential is very small at the end of inflation, r] enc i, 


V (Pend) = A 1 + cos(2arctan[c 2 77 end 3H ; 2<t2 ]) <C 1. 


(26) 
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Hence, the argument of cos should be tt. Thus, the integration constant has 
to be 

A 4 

c 2 » Ven d 3Hi2a2 ■ (27) 

The form (25) and the limit (27) will be used to derive the coupling 
function, f (r/), and then calculate the electromagnetic spectra in the de 
Sitter model approximation, in section 3. 


2.2 NI under a power law expansion 


To have a more optimal slow roll analysis that has a smooth exit from in¬ 
flation, the Hubble parameter can be written as a function of 0 , 77(0). If 
the field falls below a certain value, it starts to oscillate and then converts 
to particles in the reheating era, right after inflation. Plugging ( 6 ) and (21) 
into ( 11 ), yields 

1 _1 ( sin(OT) \ 2 

a(rj) — Io\t)\ 2c 2 v i+cos( ro ) J ^ ( 28 ) 


where, zu= arc cos 


1 



exp(—iV/C 2 


The Hubble parameter is then 


H(v) 


a'{rj) 

a 2 {rf) 


(-ii i / sin M 

' 2£ 2 Y l+COs(tZ 7 ) 

T 0 


2 



sin(-ci7) 

1+C0S('07) 


(29) 


For N > 50, and (, < 2, we have w ~ 0. Then we end up with a(r\) ~ ZoM -1 , 
and H & const, which is the same as the de Sitter model. 

However, in the more general solution, one can substitute (28) and (29) 
into (5) and solve for 0, to yield, 


0 ( 77 ) = 2 a arctan [c 2 x exp < — 


( N N \ 2 Jy_ 

e^-2C 2 + 2e^C 2 ) V ^ 2 - 2 < 2 +^ (2 c 2 A 4 

7 jv ~w \ 

6 l + ec 2 -2C 2 + 2ef 2 C 2 ) <x 


where, C 2 is the integration constant. Using the constraint, U(0(?y enc 0) -C 1, 
implies that, 


_/iv TV \ 2 N JV 

/ 0 2 (ec 2 -2( 2 + 2ec 2 C 2 J V eC ~ 2 < 2 + 2e< « 2 A 4 
6fl + eU - 2( 2 + 2eT(A a 


c 2 A> exp 


(31) 


By using (30) and (31) we can derive the coupling function, f(r/) , but it is 
very complicated and actually not needed, as the simple de Sitter approxi¬ 
mation is sufficient to investigate PMF under NI. 



3 The PMF generated in natural inflation 


Starting from Lagrangian ( 1 ), the equation of motion for the electromagnetic 
field, A^, 

[V^gg^ u g al3 f 2 (<A t) f uP ] = o, ( 32 ) 

where, g is the determinant of the flat FRW space-time metric g iw . In con¬ 
formal time, (32) can be written as 

A "i ( 77 , x) + 2 yA'j (v, x ) ~ « 2 (v) djd J Ai (77, x) = 0 ( 33 ) 

Define the function, A, (77, x) = f (77) A* (77, x), and using the quantization of 
Aj in terms of creation and annihilation operators, 6 A and b\ (k), as, 


A i (77, x)= - — (k) X [6a (k) A (77, k) e lk ' x + b j x (k) A* (77, k) e lk % 

** (2tt) X=1 

( 34 ) 

where, e % \ is the transverse polarization vector, and k = is the commoving 
wave number. Hence, Eq.( 33 ) can be written as, 

A" (77, k) + (k 2 -Y (77)) A (77, k) = 0, ( 35 ) 


n // 

where, Y (rj) = y. 

The magnetic and electric spectra can be calculated [7] respectively by, 


dpB_ 

dink 


1 

27r 2 


(77, /c)| 2 . 


(36) 


dfE _ P fc 3 
dlnfc 27 t 2 a 4 

Therefore, we need first to define the coupling function, / (77), in order to solve 
for the electromagnetic vector field, A^. We will assume that the relation 
between the coupling function and scale factor is of the power law form [7], 
/ ( 77 ) oc a a . Then, by combining (3) and (4) in the slow roll limit, 


MjiA) 

f 


( 37 ) 




oc exp 


a 


3 Mi 


f v M 

2 T' 0 )>) # 


PI 


Substituting ( 2 ) into (38) gives, 




= D sin 


0 ( 77 ) 


2qq-^ 
3 A(fpi 2 


2(7 


( 38 ) 


( 39 ) 
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where, D, is a coupling constant. As we adopt de Sitter or quasi-de Sitter 
inflationary expansion, the relation between a and 7 is a ~ —7, where 

f(v) oc ?7 7 . 

Since the power law expansion approaches de Sitter for relatively high 
N(> 50), as can be seen from (28) and (29), in this section we will only solve 
(35) in a simple de Sitter model of expansion. However, using (30) to find 
/ ( r 1 ) explicitly yields a very complicated Y ( 77 ), and an analytical solution of 
(35) cannot be found. 

The de Sitter approximation was used by [10] to investigate PMF. One 
can investigate PMF under the de Sitter model by substituting (25) into (39) 
and (35). Hence, 


f" 

Y(n) = j = 


A2arj 2+3a2H i 2 A 4 


81a 2 H 4 Mp 4 1 c 2 2 


V 


2A 4 

3a 2 H, 2 


(40) 


where, 

2A 4 

A = 3c 2 2 Mpi 2 (3a 2 Hi 2 - 2A 4 ) + <tV ct2 ^ 2 ( 9H 2 M P1 2 + 2«A 4 ) . (41) 


By using the limit (27) and the facts that, ( r/iVend ) At —1, a ~ Mpi, A ~ 
Mgut ~ 10 _ 3 Mpi[21] and the upper limit of Hi < 3.6 x 10 -5 Mpi at 95% CL 
[ 22 ], then “C 1. Thus, Y (rj) can be written as, 


Y(V) 


6ar] 2 A 4 (3a 2 H 2 — 2 A 4 ) 
81cr 2 i7j 4 Mpi 2 c 2 2 


Substituting (42) into (35), 


(42) 


- 4 " (>;%•) + 



6arj 2 A 4 (3a 2 H 2 — 2 A 4 ) \ 
81cj 2 i7j 4 Mp 1 2 c 2 2 ) 


A(r],k) 


0 


(43) 


Eq.(43) is a Bessel differential equation. Hence, A ( 77 , k ) can be written as 
A ( 77 , k) = (k V ) l/2 [C ± ( k ) J x (k V ) + C 2 (k ) J_ x (k V )\ (44) 

where % is given by 


X = 


27 + 


8a(3cr 2 J7j 2 A 4 —2A 8 ) 
C2 2 cr 2 Hi 4 M P i 2 


6V3 


10 


(45) 



In the long wavelength regime, krj <C 1 (outside Hubble radius), Eq.(44) can 
be written as 


*Akr]<Sil ij]i 



£>i (x) (~kv) X+1/2 + (x) (~kr,) I,2 - X 


(46) 


The constants, D L (y) and D 2 (x), can be fixed by using the normalization 
of A (r), k ) and other limit, Akroi (Vi k) e ~ kri / V%k- They are written as 


Di (x) = 


7 r 


3 -i 7 T( X +l/ 2)/2 


2 X+1 r (y + 1 ) cos ( 7 r(x + 1 / 2 )) 


,D 2 (x) = ^z 


7 r 


3 -i7r(x+3/2)/2 


2 1 x r (1 - x) COS (-(\ + 1/2))' 
(47) 

The magnetic spectra can be obtained by substituting of (46) into (36). 

It can be written [7] as, 


dpB = ^~(w) 4 / k \ 4+2n 

dink 27r 2 \ai7 / 


(48) 


where, 7 = x + 1/2, then n = 7 if 7 < 1/2 and n = 1 — 7 for 7 > 1/2. The 
function J-(n) can be written as, 


7-(ri) = -. 

2 2n+ ir 2 (n + 1 ) cos 2 ( 7 m) 

Similarly, the electric field spectrum can be written as, 

df>E _ g(m) 4 / k \ 4+2m 

din* 2- J l «// J 


(49) 


(50) 


where, m = 7 + lif 7 < —1/2 and m = —7 for 7 > —1/2. The function 
Q(m) can be written as, 


G(m) 


7 r 

2 3 m +3r2( m + |) cos 2 (7rm) 


(51) 


The scale invariant PMF can be achieved if the magnetic spectrum, = 
Constant. Hence, from Eq.(48), the values of 7 = {—2, 3}. The first value is 
more acceptable in generating the PMF without huge amount of backreaction 
and without assuming a small scale of inflation ( as in the case of 7 = 3) . 
Also, for a = —3, the reheating period will be to long, such that the reheating 
temperature falls to few MeV. For a = 2, the problem of strong coupling 
exists. One way to solve this problem, is to assume that, the initial coupling 
function is much less than the coupling function at the end of inflation, 
/ ( 770 ) 7C / i'i)end) ~ 1- This assumption in turns, will create a weak coupling 
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between the gauge field and charges at the end of inflation [7]. For these 
reasons, the solution, a = 2 , is adopted in this paper. 

Since, 3 a 2 A 4 2 A 8 , then % = | ^9 + - • By using the limit 

(27), we have % ~ 1/2. This value is corresponding to the value of 7 = 0 
in [7]. Hence, a scale invariant PMF cannot be generated in NI under the 
simple de Sitter model of inflation if we impose the limit (27). Calculating 
the electromagnetic spectra shows that, they are almost of the same order 
of magnitude, at krj <C 1, see Fig.l. As, 7 = 0 implies that both n,m = 
0. Hence, both magnetic and electric spectra obtained by (48-50) will be 
proportional to A: 4 , as can be seen in Fig.l. 


( d p B d p E t r 

Plnf ’ 'dhk’dfiik 1 *^ (r=0) ’ de sitterkl 1 


10 


-14 



k [Mpc 


Figure 1 : The electromagnetic spectra and the inflationary energy density, 
Pi n f, in NI at % = 1/2, rj= —20, a ~ Mpi = 1, , A = 10 -3 , Hi = 3.6 x 10 -5 ,« = 
2 and D = 1. The spectrum of electric field is of the same order of magnitude 
as the spectrum of the magnetic field for krj <C 1. At relatively high k, they 
start diverging from each other. However, the energy density of inflation 
generated by NI, pi n { is much larger than the electromagnetic energy density. 

On the other hand, if Eq.(27) limit is relaxed and the scale invariance 
condition is enforced, y = 5/2 (7 = +3, —2), then C 2 becomes 


c 2 


1 «A 4 

27 


(52) 
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Pint', 


d Pb d p E 
dink’ dink 


^X~ \ (y=-2,3), deSitterP?! 1 



Figure 2: The electromagnetic spectra and the inflationary energy density, 
Pi n f, in NI model at % = 5/2, q = —20, a ~ Mpi = 1, A = 10 -3 , Hi = 
3.6 x 10 -5 , a = 2 and D = 1. The spectrum of electric field is much greater 
than the spectrum of magnetic field for, kq <C 1. For /c min > 8.0 x 10 _7 Mpc _1 
the electromagnetic energy density can go below that of inflation. Hence, the 
backreaction problem might be avoided for k min < k -C 1. 


The coupling function can be written as, 


f(v) 


= D sin 


aretariL 


1 ctA 4 


3 H i Z f 


2 f 2 ‘ 


27 HfMy 


pi 


3 M pl 2 


(53) 


The electric and magnetic spectra in this case can be seen in Fig.2. We 
can see that at very long wavelength (k/q -C 1) the electric field spectrum 
far exceeds that of the magnetic field and the energy density of the inflation 
which is generated by NI, p In f. It may cause the backreaction problem. 
However, for P min > 8.0 x 10 _7 Mpc _1 the electromagnetic energy can go 
below that of inflation. Most of the observable scale is above k min . That 
range of k includes most of the observable scales according to Planck, 2015. 
For example, it includes the standard pivot scale, P* = 0.05Mpc _1 . Further, 
it includes some of the cut-off scale, ln(P c /Mpc _1 ) € [-12,-3], chosen by 
Planck, 2015 [22], Therefore, the backreaction problem might be avoided in 
generating PMF by the f 2 FF model in NI, under de Sitter expansion for 
1 k k m ; n . 

Likewise, plotting the electromagnetic spectra as a function of the Hubble 
parameter Hi , shows that the electric field is always greater than the magnetic 
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d p B dp £ 

dink’ dink 


^X = \ (y=-2,3),deSitter, £« I I^pi 1 



1(T 5 1(T 4 0.001 0.01 0.1 


H,[M n \ 


1 


Figure 3: The electromagnetic spectra and the inflationary energy density, 
Pi n f in NI, as a function of Hi at x = 5/2, f] = —20, cr « Mpi = 1, A = 10 -3 , 
a = 2, and D = 1. The spectrum of electric field is always much greater 
than the spectrum of magnetic field and can exceed the energy of inflation 
for H min > 1.25 x 10 -3 Mpi. This value is well above the upper limit of the 
Hubble parameter, obtained by Planck, 2015, Hi < 3.6 x 10 5 Mpi. Below 
H m i n , all electromagnetic spectra will be less than p\ n f. Hence, that might 
avoid the backreaction problem. 

field and can exceed the energy of inflation for H m ; n > 1.25 x 10 _3 Mpi.(See 
Fig.3.). This value is well above the upper limit of the Hubble parameter, 
obtain by Planck, 2015, Hi < 3.6 x 10 -3 Mpi [22], Hence, this model can be 
free from the backreaction problem. 

Similarly, plotting energy density of inflation and the electromagnetic 
spectra as a function of ( clearly shows that the backreaction problem can 
be avoided for the possible values (See Fig.4.). The value ( play an effective 
role in the characteristics of the natural inflation and its implications. For 
example, in the case of ( ;§> 1, the natural inflation behaves like quadratic 
inflation [19]. 

Finally, one can analyze the shape of the electromagnetic spectra as a 
function of A. As seen in Fig.5, there is a narrow range of A (~ 0.00874Mpi), 
at which the electric fields can even fall below the magnetic field. 

In order to decide the range of k for which the electric field energy is 
less than the magnetic fields, one can plot the electromagnetic spectra as a 
function of k, as in Fig.6. The range is k > 2.53 x 10 _3 Mpc _1 , for kr) <C 1, 
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Figure 4: The electromagnetic spectra and the inflationary energy density, 
Pi n f, in NI model as a function of ( = a/Mp\ at x — 5/2, r/ = —20, A = 10 -3 , 
a = 2, Hi = 3.6 x 10 5 Mpi and D = 1. The spectrum of electric field is 
always much greater than the spectrum of magnetic field for k 1. However, 
both electric and mgnetic spectra are much less than pi n f. The value ( 
play an effective role in the characteristics of the natural inflation and its 
implications. 
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Figure 5: The electromagnetic spectra and the inflationary energy density, 
Pi n f, in NI, as a function of A at y = 5/2, rj = —20, a ~ Mp\ = 1, a = 2, 
Z) = 1, and k = 10 _3 Mpc _1 . The spectrum of electric field falls below the 
spectrum of magnetic field around A = 0.00874Mpi. However, both of them 
are much less than p\ n f, which may avoid the backreaction problem. 

around k m j n ~ 0.0173Mpc _1 . As we choose M P i = 1( « 1 x 10 19 GeV), the 
appropriate values of A is in the order of Mqut ~ 10 16 GeV that fits with the 
results of Ref. [19]. 

On the other hand, the range of k includes most of the observable scales 
according to Planck, 2015. For example, it includes the standard pivot scale, 
k * = 0.05Mpc _1 . Further, it includes some of the cut-off scale, ln(Z c /Mpc -1 ) G 
[—12, —3], chosen by Planck, 2015 [22]. However, the relatively narrow range 
of k at which this situation is valid, may cause serious challenge to this model. 
This is so because after sufficient number of e-foldings the wave number may 
go below k < 10 _3 Mpc _1 . 

4 Summary and conclusions 

PMF can be generated by the simple inflation model / 2 FF, Eq.(l), in the 
standard models of inflation, and requires the breaking of the conformal 
symmetry of the electromagnetism. In this paper, we used the same method 
used in [7] to investigate the PMF in natural inflation. We first presented 
the slow roll analysis of the NI. Unlike the PMF on large field inflation (LFI) 
[23] and similar to PMF on R 2 -inflation [24], for sufficiently large number of 
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Figure 6: The electromagnetic spectra and the inflationary energy density, 
Pi n f, in NI as a function of k at x = 5/2, rj = —20, a ~ Mpi = 1, a = 2, 
D = 1, and A = 0.00874Mpi. The spectrum of electric field falls below the 
spectrum of magnetic field on the range of, k > 2.53 x 10 _3 Mpc _1 , at which 
the backreaction problem can be avoided. 


e-foldings, the power law inflation can lead to the same results as the simple 
de Sitter model of expansion in NI. 

We find that PMFs can in principle be generated in the NI model for all 
values of (, = a/Mp\. Under U(0) ~ 0 model of inflation, the scale invariant 
PMF is unlikely to be generated in the context of NI. That is similar to 
the case in the context of LFI [23]. However, if this constraint is relaxed, a 
scale invariant PMF can be achieved in NI. In this case, the magnitude of 
the PMF spectrum, at krj <C 1, is much smaller than the spectrum of the 
associated electric field. Changing the values of rj ,Hi, a, and D does 
not change this relation. In comparison with the inflationary energy density, 
Pi n f, in NI and the upper bound of the energy density of inflation derived 
from WMAP7, (p en d)cMB < 2.789 x 10 -lo Mpj [27], the energy of the electric 
field may exceed the energy scale of inflation at k < 8.0 x 10 _7 Mpc _1 and 
Hi > 1.25 x 10 _3 Mpi. That may prevent inflation from occurring at all. 
This is the problem of backreaction. One can conclude that for small enough 
value of k, this problem cannot be avoided in the f 2 FF model under natural 
inflation. 

On the other hand for k > 8.0 x 10 _7 Mpc _1 and Hi < 1.25 x 10 _3 Mpi 
, both electric and magnetic energy densities can fall below the inflationary 
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energy density. In this case, one can consider these values as, respectively, a 
lower bound of k and an upper bound of Hi for a backreaction-free model of 
PMF. Moreover, these scales include most of the observable ranges of k and 
Hi. 

Furthermore, there is a range of A min (~ 0.00874M P i), and k > 2.53 x 
10 _3 Mpc _1 , at which the energy density of the electric held can even fall 
below the energy density of the magnetic held. Again these values he on the 
observable range of k and the anticipated scale of Lambda. Therefore, the 
problem of backreaction can be avoided in these ranges of values. However, 
the relatively short range of k, presents a serious challenge to the viability 
of this model. One way to extend this research is to include the effect of 
reheating era and then to calculate the present value of PMF generated in 
NI as we do in the context of /A-inhation [24]. 
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